Von Kärmän vortex street in a Bose-Einstein condensate 
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Vortex shedding from an obstacle potent ial moving in a Bose-Einstein condensate is investigated. 
Long-lived alternately aligned vortex pairs are found to form in the wake, as for the von Kärmän 
vortex street in classical viscous fluids. Various patterns of vortex shedding are systematically 
studied and the drag force on the obstacle is calculated. It is shown that the phenomenon can be 
observed in a trapped System. 
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The formation of a train of alternate vortices in the 
wake past an obstacle, known as the von Kärmän vor- 
tex street, is a ubiquitous and intriguing phenomenon 
in fluids. Since the pioneering experimental study 
by Benard [H and theoretical consideration by von 
Kärmän [2], numerous studies have been made on the 
phenomena of vortex street formation jsf. 

The behavior of a viscous fluid flowing past an ob- 
stacle is determined by the Reynolds number Re, which 
is a dimensionless parameter that includes the kine- 
matic viscosity. The wake of a cylinder is steady for 
Re below around 50, and a vortex street emerges for 
10^ ^ Re ^ 10^, which becomes turbulence for a larger 
Re. This implies that the viscosity plays an Import ant 
role in vortex street formation in classical fluids. For su- 
perfluids, however, the Reynolds number cannot be de- 
fined because of the absence of viscosity. Moreover, the 
vortex quantization makes superfluid dynamics quite dif- 
ferent from classical fluid dynamics. Therefore, it is not 
obvious whether instability of the wake and subsequent 
vortex street generation occur in superfluids. According 
to von Kärmän's theory 0, a vortex street is expected 
to be very long-lived in inviscid fluids, once it is created. 

In this Letter, by numerically solving the Gross- 
Pitaevskii (GP) equation, we show that long-lived al- 
ternately aligned vortex pairs are formed in the wake 
of an obstacle potential moving in a Bose-Einstein con- 
densate (BEC), which is similar to the von Kärmän vor- 
tex street in classical fluids. Mean-field analysis of Sys- 
tems of a BEC with a moving potential has been per- 
formed by many authors from the viewpoints of drag 
force [Ü, , vortex dynamics near the cylinder Q , critical 
velocity [Sil, supersonic flows [9|,110], and multicomponent 



the GP equation given by 



Systems [ll|, ll2|. However, vortex street formation was 



not found in these studies, probably because the param- 
eter region for which a vortex street emerges is narrow, 
as shown later. 

We consider a BEC of atoms with mass m and an 
obstacle potential V moving in the —x direction at a 
velocity v. In the mean-field theory, the condensate is 
described by the macroscopic wave function obeying 
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where g = AnTi^a/m with a being the 5- wave scattering 
length of the atoms. We employ a Gaussian potential 
with peak strengt h and radius d moving in the —x 
direction at a velocity v d^s V = Vbexp{ — [(x + vt^ + 
y'^]/d?}. Normalizing space and time by h/{mgnoy^'^ 
and fi/{gno)^ where uq is the atom density far from the 
potential, we can eliminate the interaction parameter g 
from Eq. ([1]). We numerically solve Eq. ^ in two di- 
mensions under the periodic boundary condition using 
the pseudo-spectral method. The initial state is the sta- 
tionary state of Eq. ([T]) with v = plus a small amount 
of noise to break the symmetry. 

Figure [T] shows typical wakes fiowing past an obsta- 
cle potential with Vo/{gno) = 100. For a sufficiently 
small velocity the fiow around the obstacle is a steady 
laminar flow and no quantized vortex is created. When 
the velocity v exceeds a critical velocity, which depends 
on the interaction strength and the shape of the poten- 
tial, vortex- antivortex pairs are created [5]. The crit- 
ical velocity for the vortex creation is of the order of 
the speed of sound (^no/m)^/^. When a created vortex- 
antivortex pair separates from the potential, the flow 
velocity around the potential again exceeds the critical 
velocity and a subsequent vortex- antivortex pair is cre- 
ated. A train of vortex- antivortex pairs is thus generated 
behind the potential. Since a Symmetrie double row of 
vortices is unstable [4], the vortex pairs are dislocated 
sinuously as shown in Fig. [T](a). Nore et al. [l3| showed 
that such a staggered vortex pattern is formed if a double 
row of vortices is prepared with an appropriate pertur- 
bation. Since a pair of point vortices with circulations 
±h/m {h: Planck's constant) and distance d moves in 
the direction perpendicular to a line between the pair 
at a velocity h/{md) [4], the alternately inclined vortex 
pairs move in two directions [white arrows in Fig. [T] (a)], 
forming a V-shaped wake as in Fig. [T] (a). The diver- 
gence of the wake is signiflcant for large which forms 
a pattern similar to supersonic flow [9|, llO| . 

Figure [T] (b) shows the main result of this study. The 
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(a) V = 2.4, d / ^ = 0.04 (b) v = 2.6, d / ^ = 0.05 (c) v = 3.0, d / ^ = 0.05 




FIG. 1: (Color) Density and phase distributions of a condensate past an obstacle potential. The velocity and potential width are 
(v, d/0 = (2.4, 0.04) in (a), (2.6, 0.05) in (b), and (3.0, 0.05) in (c), where v = v[10^m/{gno)]^^^ /{27i) and ^ = ;i[10V(^^^o)]^/^ 
The white arrows in (a) indicate the directions in which the vortex-antivortex pairs move. The density is normahzed by uq. 
The field of view is 6^ x 3^. In the numerical calculation, a 32^ x 8^ space is discretized into 4096 x 1024. 



significant difference from Fig. [T] (a) is that the vortices 
in a pair created by the obstacle potential at a time have 
the same circulation. Since two point vortices having 
the same circulation h/m rotate around their center at 
an angular frequency 2?i/{mdP) without changing their 
distance [4], the created vortex pairs in Fig. [T] (b) re- 
main bound and rotate. The pairs with opposite cir- 
culations are alternately released from the obstacle po- 
tential to form a train of vortex pairs, resembling a von 
Kärmän vortex street. In contrast to the vortex arrange- 
ment originally considered by von Kärmän, in which iso- 
lated point vortices are aligned, the vortex pairs consti- 
tute the vortex street in the present case. We find from 
Fig. [1] (b) that the distance between the two vortex rows 
\s h 0.24^ and the distance between two pairs in a row 
is :^ 0.87^ on average, and hence 6/^ :^ 0.28, where ^ = 
/i[10^/(m^no)]^/^. This ratio is in good agreement with 
the stability condition of von Kärmän's vortex arrange- 
ment hß = tt'^ cosh"^ 0.28 0,13. In fact, the vor- 

tex street in Fig. [T] (b) survives at least t > 10^h/{gno). 
The vortex street in Fig. [T] (b) moves in the —x direction 
at a velocity O.U{gno/m)^^'^ O.Sh/{V2im). The 
velocity of von Kärmän's point vortices, in which each 
vortex has a circulation 2/i/m, is h/{V2£m) 0, Ü]. For 
large d and the periodicity in the wake seems to dis- 
appear [Fig. [T] (c)]. We have numerically confirmed that 
similar wakes are also obtained for a disk-shaped poten- 
tial {V = oo foT r < d and V = for r > d). 

Figure[2]shows the dynamics of vortex street formation 
just behind the obstacle potential. The pairs of vortices 
are released obliquely backward left and right with alter- 
nate circulations. This vortex street formation behavior 
is quite different from that in classical viscous fiuids, in 
which a pair of eddies in the wake becomes unstable and 
grow into a vortex street in the downstream region. In 
contrast, we find from Fig. [2] that the alternate vortices 
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FIG. 2: (Color) Serial snapshots of the density profiles for 
the Parameters in Fig. [1] (b) in the frame moving with the 
potential. The time interval is 10h/ (gno). The field of view 
is 2^ X <f . The color scale is the same as that in Fig. [Ü 



are directly created at the obstacle, implying that the 
mechanism of vortex street formation in Fig. [2] may be 
different from that in classical fiuids. 

We systematically performed numerical simulations for 
various values of d and v to determine the parameter 
regions for the types of wakes in Fig. [TJ Figure [3] shows 
a rough sketch of each parameter region. The regions 
of the periodic behaviors shown in Figs. [T] (a) and [T] (b) 
are located between the regions of steady laminar fiow 
(white region in Fig. [3]) and irregulär fiow (red). We note 
that the parameter region for vortex street formation is 
rather restricted, 0.04 < d/^ < 0.13 and 1.9 < v < 2.8, 
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FIG. 3: (Color) Dependence of the patterns of wakes on the 
normalized Gaussian width d/^ and velocity v of the poten- 
tial. The green, blue, and red regions correspond to the flow 
patterns shown in Figs. [U (a), [1] (b), and[T] (c), respectively. 
The white region corresponds to stationary laminar flow. 
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FIG. 4: (Color online) Time evolution of the normalized drag 
force F = Fnm~^^\lO~^ gno)~^^^ . The solid and dashed 
lines show Fx and Fy. The parameters used in (a)-(c) are 
the same as those in Figs. [T] (a)-(c), respectively. Time is 
normalized as t = 10~^tgno/fi. 



where v = v[10^m/ {gno)]^^'^ / {27r). This is in contrast 
with classical fluids, in which the von Kar man vortex 
Street emerges for a wide ränge of Reynolds number. 

Figure |4] shows the drag force on the obstacle Poten- 
tial given by F = / drip* The initial state is 
the stationary state for v = 0. At t = the potential 
Starts to move at a velocity v. Figure |4](a) corresponds 
to the vortex-antivortex pair creation in Fig. [1] (a). For 
t ^ 20, Fx oscillates while 0, indicating that the 
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FIG. 5: (Color) Density proflle of a BEC of ^'^Rb atoms 
conflned in a harmonic potential at t = 330 ms. The Gaussian 
Potential is initially located at the marked position (x) and 
Starts to move in the —x direction at a constant velocity. The 
density and phase proflies in the dashed region are magnifled 
in the lower panels. The fleld of view is 430 x 43 /xm for the 
main panel and 86 x 21 /im for the lower panels. The color 
bar is normalized by Nmujy/Ti. 



vortex-antivortex pairs are shed from the potential sym- 
metrically. For t ^ 20, the vortex pairs begin to incline as 
in Fig. [T](a), and Fy also Starts to oscillate. Figure |4](b) 
corresponds to the vortex street formation, where both 
Fx and Fy oscillate for t ^ 20. The oscillation in Fy in 
Figs. m (a) and |4] (b) is due to the alternate shedding of 
vortices and hence its frequency is half the vortex shed- 
ding frequency, i.e., the frequency of F^. It is interesting 
to note that Fy in Fig.|4](c) oscillates, even though Fig.[T] 
(c) does not seem to have periodicity. A similar phe- 
nomenon is also observed in classical fluids with a large 
Reynolds number [13]. 

Next, we study a realistic System conflned in a trap- 
ping potential. A BEC stirred by a moving obstacle po- 
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tential has been studied experimentally in Refs. 
We consider a Situation in which a BEC of ^^Rb atom is 
conflned in a harmonic potential m(u;^x^-hco'^^^+cc;^z^)/2 
with {üJx,üJy,üJz) = 27r x (2.5,25, 1250) Hz. The number 
of atoms is A/" = 10^ and the s-wave scattering length 
is lOOttß with ttß being the Bohr radius. The conden- 
sate is tightly conflned in the z direction and the System 
is effectively two dimensional. The obstacle potential is 
produced by a blue-detuned Gaussian laser beam with 
d = 2.16 /im and Vq = lOOficj^y, which is initially located 
at X = 130 /im and y = at t = and moves in the —x 
direction at a velocity v = 0.68 mm/s for t > 0. Fig- 
ure [5] shows the density and phase proflies at t = 330 
ms. We flnd that a vortex street is gener ated behind the 
moving potential, as in Fig. [T] (b), which conflrms that 
our flnding can be experimentally observed in an inho- 
mogeneous System. The ratio hß in Fig.[5]is about 0.22, 
which deviates from 0.28 due to the flnite size effect. 

A future prospect of this study is to provide a mathe- 
matical description of the symmetry breaking instability 
and clarify the mechanism of the alternate creation of the 
vortex pairs as shown in Fig.[2l In classical viscous fluids, 
the transition from stationary flow to a vortex street is 
characterized by the Hopf bifurcation and described by 
the Stuart-Landau model [18|]. A similar approach may 
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be applied to the present phenomenon. Various shapes of 
obstacle potentials and three dimensional dynamics also 
merit further study. 

In conclusion, we have shown that pairs of quantized 
vortices shed from an obstacle potential moving in a BEC 
alternately align and survive for a long time [Figs. [T] (b) 
and [2] , which closely resembles the von Kärmän vortex 
Street in classical fluids. We have obtained the parame- 
ter region in which a vortex street emerges (Fig. [3]) and 
calculated the drag force (Fig. HJ. We have shown that 
vortex street formation can be observed in a trapped 
BEC disturbed by a blue-detuned laser beam (Fig. [5]). 
Since the von Karmän vortex street typifies the great di- 
ver sity of classical fluid dynamics, the emergence of the 
von Kärmän vortex street in a BEC implies that a rieh 
variety of phenomena is still unrevealed in quantum hy- 
dro dynamics. 

We thank T. Miyazaki and N. Takahashi for their valu- 
able comments. This work was supported by the Min- 
istry of Education, Culture, Sports, Science and Tech- 
nology of Japan (Grants-in-Aid for Scientific Research, 
No. 17071005 and No. 20540388). 



[1] H. Benard, C. R. Acad. Sei. Paris 147, 839 (1908); 147, 
970 (1908). 

[2] T. von Kärmän, Naehr. Ges. Wiss. Göttingen, Math. 

Phys. Kl. 509 (1911); 547 (1912). 
[3] For review, see for example, C. H. K. Williamson, Annu. 

Rev. Fluid Mech. 28, 477 (1996). 
[4] For example, H. Lamb, Hydrodynamics, 6th ed. (Dover, 



New York, 1945), Chap. VII. 
[5] T. Frisch, Y. Pomeau, and S. Rica, Phys. Rev. Lett. 69, 
1644 (1992). 

[6] T. Winiecki, J. F. McCann, and C. S. Adams, Phys. Rev. 

Lett. 82, 5186 (1999). 
[7] C. Nore, C. Huepe, and M. E. Brächet, Phys. Rev. Lett. 

84, 2191 (2000). 
[8] J. S. Stießberger and W. Zwerger, Phys. Rev. A 62, 

061601(R) (2000). 
[9] G. A. El, A. Gammal, and A. M. Kamchatnov, Phys. 

Rev. Lett. 97, 180405 (2006). 
[10] I. Carusotto, S. X. Hu, L. A. Collins, and A. Smerzi, 

Phys. Rev. Lett. 97, 260403 (2006). 
[11] H. Susanto, P. G. Kevrekidis, R. Carretero-Gonzälez, 

B. A. Malomed, D. J. Frantzeskakis, and A. R. Bishop, 

Phys. Rev. A 75, 055601 (2007). 
[12] A. S. Rodrigues, P. G. Kevrekidis, R. Carretero- 
Gonzälez, D. J. Frantzeskakis, P. Schmeicher, T. J. 

Alexander, and Yu. S. Kivshar, Phys. Rev. A 79, 043603 

(2009). 

[13] C. Nore, M. E. Brächet, and S. Fauve, Physica D 65, 154 
(1993). 

[14] A. Roshko, J. Fluid Mech. 10, 345 (1961). 

[15] C. Raman, M. Köhl, R. Onofrio, D. S. Durfee, C. E. 

Kuklewicz, Z. Hadzibabic, and W. Ketterle, Phys. Rev. 

Lett. 83, 2502 (1999). 
[16] R. Onofrio, C. Raman, J. M. Vogels, J. R. Abo-Shaeer, 

A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. 85, 

2228 (2000). 

[17] T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis, 

and B. P. Anderson, arXiv:0912.3773 
[18] C. Mathis, M. Provansal, and L. Boyer, J. Phys. Paris 

Lett. 45, 483 (1984); M. Provansal, C. Mathis, and L. 

Boyer, J. Fluid Mech. 182, 1 (1987). 



